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ABSTRACT 

We conduct a detailed comparison of broad-band spectral energy distributions of six 
z > 5.5 galaxies against galaxies drawn from cosmological hydrodynamic simulations. 
. We employ a new tool called Spoc, which constrains the physical properties of ob- 

served galaxies through a Bayesian likelihood comparison with model galaxies. We 
first show that SPOC self-consistently recovers the physical properties of a test sample 
of high-rcdshift galaxies drawn from our simulations, although dust extinction can 
yield systematic uncertainties at the w 50% level. We then use Spoc to test whether 
, our simulations can reproduce the observed photometry of six z > 5.5 galaxies drawn 

t~^. ■ from the literature. We compare physical properties derived from simulated star forma- 

C^) ' tion histories (SFHs) versus assuming simple models such as constant, exponentially- 

^sO . decaying, and constantly rising. For five objects, our simulated galaxies match the 

observations at least as well as simple SFH models, with similar favored values ob- 
(-h \ tained for the intrinsic physical parameters such as stellar mass and star formation 

O rate, but with substantially smaller uncertainties. Our results are broadly insensitive 

to simulation choices for galactic outflows and dust reddening. Hence the existence 
of early galaxies as observed is broadly consistent with current hierarchical structure 
formation models. However, one of the six objects has photometry that is best fit by 
a bursty SFH unlike anything produced in our simulations, driven primarily by a high 
iiT-band flux. These findings illustrate how SPOC provides a robust tool for optimally 
utilizing hydrodynamic simulations (or any model that predicts galaxy SFHs) to con- 
strain the physical properties of individual galaxies having only photometric data, as 
well as identify objects that challenge current models. 



o 

to 

in 
c3 



Key words: galaxies: formation, galaxies: evolution, galaxies: high-redshift, cosmol- 
ogy: theory, methods: numerical 



1 INTRODUCTION 

Over the last few years, observations of galaxies at 
z ~ 6 have opened up a new window into the reion- 
ization epoch, turning it into the latest frontier both 
for observational an d theoretical studies o f galax y for- 
mation. Planned (|Gonzalez-Serrano et al.l 1200 5) and 
existing wide-area narrowband searches for z > 5 . 5 ob- 
jects such as th e Sub aru Deep Field (jAiiki et alj 120061 : 
20061). the Large Area Lyman-Alpha 
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targeted searches near lensing caustics in galaxy clus- 



ters dKneib et al] |2004 iHu et al] I2OO2T ) and occasionally 
serendipity (IStern et al.l l2005f i to uncover star-forming 
galax ies from the reioni zation epoch in significant numbers 
(see iBerger et all [2006, for a listing of spectroscopically- 
confirmed z > 5 galaxies). 

A question immediately raised by this new stream of ob- 
servations is, what are the physical properties of these early 
galaxies? Optimally, one would determine properties such as 
the stellar mass, star formation rate, and metallicity directly 
from high-quality spectra, but at present this is infeasible for 
such faint systems. Hence properties must be inferred from 
photometry alone, occasionally augmented by emission line 
information. This requires making some poorly constrained 
choices for the intrinsic galaxy properties. A commonly ap- 
plied method known as Spectral Energy Distribution (SED) 
fitting involves generating an ensemble of population syn- 
thesis models under a range of assumptions for the intrinsic 
nature of the object, and then finding the set of assumptions 
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that be st reproduces a Riven galaxy's obser ved photome- 
try (e.g. lBemtej|2000l ; lKauffmann et al.ll2003h . The physical 
properties that yield the lowest \ 2 model are then forwarded 
as the most probable values, sometimes wi th little attention 
to sta tistical uniqueness or robustness fsee lSchaerer fc Pelldl 
2005, for a nice exploration of such issues). 

Amongst the various assumptions used in SED fitting, 
the one that is often least well specified and produces the 
widest range in final answers is the galaxy's star formation 
history (SFH). With no prior information, common practice 
is to use simple SFHs with one free parameter such as con- 
stant, single-burst, or exponentially-decaying, which in ag- 
gregate are assumed to span the range of possible SFHs for a 
given galaxy. Indeed, in most cases all one-parameter SFHs 
yield plausible results, though the parameters obtained and 
quality of fits in each case can vary significantly. If it were 
possible to narrow the allowed range of SFHs through inde- 
pendent considerations, physical parameters could in prin- 
ciple be more precisely determined. 

One approach for constraining SFHs a priori is to incor- 
porate information from currently favored hierarchical struc- 
ture formation models. As we will discuss in this paper, hy- 
drodynamic simulations tend to produce a relatively narrow 
range of star formation histories for early galaxies. Their 
galaxies' SFHs tend to follow a generic form at these early 
epochs, best characterized as a constantly-rising SFH. This 
form is broadly independent of cosmology, feedback assump- 
tions, or other ancillary factors, and is furthermore distinct 
from any one-parameter models commonly used today. A 
primary aim of this paper is to test whether this relatively 
generic SFH form is consistent with observations, and if so, 
what the implication are for the physical properties of high- 
redshift galaxies. 

Despite impressive recent successes in understanding 
cosmology and large-scale structure in our Univer se (e.g. 
ISpergel et all [2006; Spr ingel. Frenk. fc White! 12006), many 
uncertainties remain in our understanding of galaxy forma- 
tion. Several recent papers have tested models of high-z 
galaxy formation by comparing them to observed bulk prop- 
erties such as luminosity functions at rest-frame UV and Lya 
wavelengths. These comparisons have shown that such mod- 
els are broadly successful at reproducing observations, under 
reasonable assumptions f or poorly constrained parameters 
such as dust extinction (Somerville et al. 200 ll; Idzi et al 



20041; iNight et al.l 120051 ; iFinlator et all l200rj ; iDave et al 



2006al ). While this broad success is encouraging, it is sub- 



ject to some ambiguousness in interpretation, because the 
properties of individual galaxies are not being compared in 
detail. One could envision situations in which a model repro- 
duces an ensemble property of galaxies but not the detailed 
spectra of individual objects. As an example, it was for- 
warded by iKolatt et all (f999) that Lyman break galaxies 
at z ~ 3 are actually merger-driven starbursts, in contrast 
to many other models predicting them to be large quiescent 
objects. Despite quite different SFHs, both models repro- 
duced many of the same bulk properties such as number 
densities and clustering statistics. For z > 6 galaxies where 
statistics are currently poor, such degeneracies can hamper 
interpretations of bulk comparisons of observations to mod- 
els. 

A complementary set of constraints on galaxy forma- 
tion models may be obtained by comparing models to the 



individual spectra of observed galaxies. In practice, for high- 
z systems, photometry over a reasonably wide set of bands 
must substitute for detailed spectra. Such comparisons of 
models to data would move towards more precise and sta- 
tistically robust analyses that do not rely on having a large 
ensemble of objects. This last aspect is critical, because the 
very earliest observed objects that may provide the greatest 
constraints on models will in practice always be few in num- 
ber and detected only at the limits of current technology. 

In short, what is desireable would be a tool to com- 
pare models and observations of high-redshift galaxies that 
(1) employs reasonably generic predictions of current galaxy 
formation models; (2) provides a quantitative and robust 
statistical assessment of how well such models reproduce 
observations; (3) yields information on the physical proper- 
ties of galaxies under various assumptions; (4) obtains such 
information based solely on observed photometry; and (5) 
does all this on a galaxy-by-galaxy basis rather than relying 
on having a large statistical sample of observed galaxies. 

In this paper we introduce such a tool, called 
Spoc (Simulated Photometry-derived Observational Con- 
straints). Spoc takes as its input the photometry (with 
errors) of a single observed galaxy along with an ensem- 
ble of model spectra drawn either from simulations or gen- 
erated using one-parameter SFHs. The output is proba- 
bility distributions of physical parameters derived using a 
Bayesian formalism, along with goodness-of-fit measures for 
any given model. The probability distributions give quanti- 
tative constraints on the galaxy's physical properties, while 
the goodness-of-fit can be used to discriminate between 
models and determine whether a given model (be it sim- 
ulated or one-parameter SFHs) is able to provide an accept- 
able fit to that galaxy's photometry. 

After introducing and testing Spoc, we apply it to a 
sample of six z > 5.5 galaxies from the literature that have 
published near-infrared photometry. We show that in five 
of six cases, the simulated galaxies fit observations at least 
as well as one-parameter SFHs. Since there is no guarantee 
that simulations produce galaxy SFHs that actually occur 
in nature, the fact that good fits are possible shows that 
the existence of the majority of observed z > 5.5 galaxies is 
straightforwardly accommodated in current galaxy forma- 
tion models. However, in one case, we find that simulated 
galaxies provides a much poorer fit than can be obtained 
with one-parameter SFHs, as burstier SFHs provide a much 
better fit than can be obtained from any simulated galax- 
ies. At face value, this implies that our simulations cannot 
yet accommodate the full range of observed galaxies, and 
that some physical process may be missing, although we 
will explore alternate interpretations. For each galaxy we 
also present the best-fit physical parameters, with uncer- 
tainties, obtained using each model SFH. The simulations 
provide significantly tighter constraints than the full range 
of one-parameter SFHs, as expected based on their relatively 
small range of SFHs produced. These values can therefore 
be regarded as predictions of our simulations that may be 
tested against future observations. 

§ [2] introduces Spoc, detailing our Bayesian formalism 
and discussing systematic uncertainties. §[3]presents the sim- 
ulations and the one-parameter models that will be used as 
the template library for Spoc. §|4]discusses what drives the 
inferred physical properties in the context of our simulations, 
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and shows that Spoc accurately recovers the physical prop- 
erties of simulated galaxies. §[5]explores the best-fit parame- 
ters of one observed reionization-epoch galaxy in detail, and 
compares with results from traditional one-parameter SFH 
models. §[6] repeats the previous comparison for a larger set 
of observed galaxies, highlighting the variety of interesting 
results that Spoc obtains. Finally, in § [7] we present our 
conclusions. 



2 METHODOLOGY OF Spoc 
2.1 SED Fitting 

Pedagogical explanations of SED fittin g techniques have 
been presented elsewhere (|Benftej |2000| ; iKauffman n et ah! 
2003), so we refer the reader there for more detailed discus- 
sion of those aspects. Here we provide some basic insights 
and notes. 

Clearly, the amount of physical information that can 
be inferred from available data depends on the quantity and 
quality of the data. For some high-z galaxies, only narrow- 
band photometry and res t-frame ultraviolet (UV) spec- 
trosc o py are availab l e (e.g . JCubv et al. 2003; Kodaira et al.l 



2003; Rhoads et al. 2003; Kurk et al. 2004; Rhoads ct al. 



2004 ; IStern et al.ll2005l ; IWestra et al-lbOOSp T For others, 



emission-line m easurement and 1-3 rest-UV broad bands a re 
available (e.g., Nagao et alj 2004 ■ IStanwav et alj 2004a 
iNagao et al.l 120051 ; IStiavelh et all 120051 ; iHu et all 1200- 



Studies employing the Lyman dropout technique in the op- 
tical must further co ntend with the possible presence of low- 
redshi ft interlopers dDickinson et al.l |2004| ; iBouwens et al.l 
and large uncertainties from dust extinction. Nev- 
ertheless, some interesting constraints can be placed on the 
underlying ph ysical properties of the sources from sole ly 
rest-UV data fProrv et al.ll2005l ; lGwvn fc Hartwickll2005l ). 



With the addition of rest-frame optical data, e.g. from 
Spitzer's Infrared Array Camera (IRAC), it becomes pos- 
sible to obtain simultaneous constraints for the stellar 
mass, star formation rate (SFR), dust extinction, and red- 



smrt using spectral cnerg 
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2004h. The uncertainties inherent in such analyses primarily 



stem from a poor constraint on the age of the galaxy's stel- 
lar population, because the relationship between age and the 
strength of the tel ltale Balmer break depends on the form of 
the assumed SFH l|Papovich et al.ll200ll ; IShaplev et al1l2005l: 
Figure [8]). This age uncertainty propagates via a host of de- 
generacies into increased uncertainties in the inferred stel- 
lar mass, SFR, metallicity, and dust extinction, if no priors 
are assumed on these quantities. Additional uncertainties 
arise fro m the unknown form of the appropriate template 
SED (e.g. ISchaerer fc Pelloll2005h and the treatment of stel- 
lar evolution as sumed by the chosen population synthesis 
models (see e.g. iMaraston et al.l[200r3 ). Still, SED fitting of- 
fers the most promising approach for determining the phys- 
ical properties of individual high-z galaxies. 

Given this, how can one employ simulations to improve 
constraints on SED fitting? One can view a numerical sim- 
ulation as producing a Monte Carlo sampling of parameter 



space such that the frequency with which a given set of phys- 
ical parameters ought to occur is proportional to the number 
of galaxies in the simulation that are characterized by that 
set of parameters. In essence, numerically-simulated galax- 
ies provide "implicit priors" for SED fitting, i.e. solutions 
that are a priori weighted more heavily because they occur 
more frequently. 

The underlying assumption is that simulated galaxy 
SFHs represent those occuring in nature. This is by no 
means guaranteed, and indeed whether Spoc provides an 
acceptable fit to a given galaxy constitutes a stringent test 
of the simulation, because a galaxy's spectrum encodes in- 
formation about its full SFH. This is the manner in which 
Spoc can provide a test of galaxy formation models based 
on individual systems. 

2.2 The SPOC Equation 

We now summarize the Bayesian statistical method em- 
ployed in SPOC. Our goal is to constrain the stellar mass, 
SFR, mean stellar metallicity, age, dust extinction, and 
redshift (M*, M„, Z, , t, Ay, and z, respectively) based 
on available measurements D. According to Bayes' The- 
orem, the probability p that the measurements D corre- 
spond to a galaxy with the intrinsic physical parameters 
<j) = (A/*, M*, Z*, t, Ay, z) (where a hat indicates a particu- 
lar value of a parameter) is given by 



p{<f> | D)<xp{<t>)p{D | <t>). 



(1) 



The prior p(4>) indicates the relative a priori probability 
that a randomly selected galaxy has this particular com- 
bination of parameters, and the likelihood p(D \ (f>) indi- 
cates the probability of obtaining the measurements D for a 
galaxy characterized by the parameters <f>; for a given model 
galaxy and data set D this is assumed to be proportional 
to e x ' . Any information regarding the expected distribu- 
tions of physical properties of the observable galaxies (such 
as the stellar mass function) or relationships between these 
properties (such as a mass-metallicity relation) can be taken 
into account via a contribution to the prior, and will gener- 
ally give rise to more precise — and possibly more accurate — 
constraints. 

In this work, we assume uniform priors on z and Ay, 
and we do not assume any dependence between Ay and 
the other intrinsic physical properties. We introduce an ad- 
ditional prior p(sim) to account for any other priors. For 
example, when matching observed galaxies against model 
galaxies derived from the outputs of two cosmological sim- 
ulations that span different comoving volumes, p(sim) rep- 
resents the ratio of the simulation volumes. After several 
applications of the product rule, we obtain 



p(4>) oc p(M* , M t , Z, , i | z)p(sim)p(D 



(2) 



This is the fundamental equation that Spoc evaluates. 
Generically, one would use Equation [2] by beginning with 
a set of models that uniformly samples the relevant pa- 
rameter space and then guessing the form of the prior 

p(M*, M*, Z*,t | z), which now encodes the assumed distri- 
bution of intrinsic physical properties of galaxies as a func- 
tion of redshift. In the high-redshift literature, where little is 
known about the intrinsic physical properties of the galaxies, 
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it is common to neglect priors altogether (or, equivalently, 
to choose the model with the lowest x 2 ) or even to intro- 
duce them accidentally by not sampling parameter space 
uniformly. The difference between this work and that of pre- 
vious authors is that we account for this prior implicitly by 
using numerically simulated galaxies as the model set. 

To see how this works, consider how one would use equa- 
tion [2] in practice. For simplicity, suppose that we wished to 
constrain a galaxy's stellar mass and that the mass could 
only fall within one of two ranges. If we omitted priors and 
assumed that the models sample stellar mass uniformly, 
then the probability that the galaxy's mass falls within a 
given range would be given by Ae~ Xi ^ 2 , where the sum 
is taken over all models whose mass lies within that range 
and the normalization A is chosen so that the sum taken 
over all models in both ranges equals unity. If we believed 
that galaxies with masses in one range were, say, twice as 
common (and therefore a priori twice as likely to be the 
right answer) as galaxies with masses in the other range, we 
could account for this via an explicit prior by changing the 
sum to APie~ Xi/l2 where Pi — 2 for models in the more 
common range and 1 for models in the less common range 
(with A of course rescaled). It is clear that an equivalent 
method to employing this explicit prior would be to gener- 
ate twice as many models in the more common range, result- 
ing in twice the probability of selecting one of these models. 
Generalizing this idea, one can view simulated galaxies as 
a Monte Carlo sampling of parameter space that naturally 
produces more models with parameters that are more com- 
monly found. Hence by taking a set of simulated galaxies, 
generating a library by resampling this set with parameters 
having uniform priors (namely, Ay and z), and using that 
library to discretely sample the probability distribution in 
the right-hand side of equation [5] one can solve equation [2] 
effectively incorporating the implicit priors given by the sim- 
ulated galaxies. This is in essence the Spoc algorithm. 



2.3 Systematic Uncertainties in Using Simulated 
Galaxies 

A major difficulty with the Spoc approach is that there is 
no guarantee that the simulation predicts the correct distri- 
bution of intrinsic properties of galaxies; in Bayesian terms, 
the priors could be wrong. On some level this is bound to be 
the case as we do not account for every process that could in 
principle affect galaxies at this epoch; indeed, no model cur- 
rently does. However, our goal is to determine whether our 
treatment is sufficient to account for current observations. If 
not, then the failures indicate needed improvements to the 
model. If our treatment can account for current observa- 
tions, then the constraints that we derive may be regarded 
as physically-motivated predictions, subject to verification 
when more constraining data become available. 

The two greatest uncertainties for the input physics 
present in our current simulations are (1) numerical 
resolution — manifested either as an inability to account for 
physical processes that occur on scales that are too small or 
too rapid (e.g. merger-induced starburst) or as a lack of nu- 
merical convergence — and (2) the prescription for superwind 
feedback. In § 14.31 we use a simple convergence test to ar- 
gue that our results do not suffer from numerical resolution 



limitations. As to our treatment for outflows, we can esti- 
mate the extent of any resulting systematics by comparing 
results from our three different outflow simulations. While 
this does not span the full range of possible feedback mecha- 
nisms, the fact that (as we show in § 15. ip most of the best-fit 
parameters are insensitive to the choice of wind prescription 
suggests that outflows do not noticeably alter typical SFHs 
at a given stellar mass. 

On the other hand, if there are significant physical 
processes affecting galaxy SEDs that are not accounted 
for by our simulation or population synthesis models, then 
our simulated galaxies may fail to reproduce the observed 
spectra, or they may mistakenly model nonstellar contri- 
butions to the observed SED as starlight. Among the pos- 
sibilities here are active galactic nuclei (AGN), incorrectly 
modeled ther mally pulsating asym ptotic giant branch (TP- 
AGB) stars (IMaraston et al.l 120061 ). emission lines, and an 
inappropriate treatment of dust or IGM absorption. We will 
argue in fj6]that significant AGN contamination is unlikely 
for the high-redshift objects we will consider here. The con- 
tribution of TP-AGB stars is also unlikely to be important 
partly because we do not model measurements from bands 
redder than / in the rest-frame, and partly because at z ~ 6 
less than half of the existing stellar mass is more than 200 
Myr old (Table [2}. Emission lines and incorrectly-modelled 
IGM abso rption could in princip l e affect our result s at the 
10% level l|Schaerer fc Pelldll2005l ; [Egami et al.ll2005l ). These 
effects are expected to be similar for the various SFHs inves- 
tigated because we use the same population synthesis mod- 
els to model the stellar continuum in eac h case. Regarding 
dust, we have found, in agreement with ISchaerer fc Pelldl 
(2005), that our results are relatively insensitive to the form 
of the dust law that we consider (see § I4.1J1 . Thus, for the 
preliminary study in this paper we ignore all of these effects. 

Another possible problem is that galaxy classes that 
are rare in reality are likely to be rare in the simulations. 
Accordingly, if the comoving volume from which the cata- 
log of simulated comparison galaxies is drawn is sufficiently 
small that a simulated analogue to an observed rare ob- 
ject is neither expected nor found, that object cannot di- 
rectly constrain the model. For example, our simulations 
produce no galaxies massive enough at z > 6 to fit HUDF- 
JD 2, the putative 6 x 10 1:L Afo object at z ~ 6.5 reported 
by iMobasher et all l|2005l '). Alt hough this particula r object 
is likely to be at a lower redshift l|Dunlop et al.ll2006T ). it does 
illustrate limitations imposed by simulation volume, which 
could also impact const raints on rare cla sses such as sub- 
millimeter galaxies (e.g. ISmail et al.ll2004h (alternatively, if 
such objects are indeed common at z ^ 6 then they represent 
a challenge to our simulations). In principle one could work 
around this issue by running larger-volume simulations or 
by deriving the priors from the simulations and then resam- 
pling parameter space by hand. In lieu of these approaches, 
the simulations utilized must have comoving volumes com- 
parable to the effective volume of the survey in which the 
object was found. 

It may appear overly ambitious to attempt to constrain 
6 (or more) seemingly independent parameters for a galaxy 
for which fewer than 6 measurements are available. How- 
ever, cosmological simulations allow us to do this because 
they generically predict that galaxies' intrinsic physical pa- 
rameters are manifestly not independent; there are tight 
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predicted correlations between, for example, stellar mass on 
the on e hand and star formation rate and met allicity on the 
other (jFinlator et alJl200d : I5ave et alJl2006al ). 

In summary, using simulated galaxies to estimate phys- 
ical properties is only valid when the dominant emission 
mechanism is star formation, and when other uncertainties 
can be carefully analyzed and shown to be negligible. For 
galaxies at high redshift, such as the ones we consider in 
this paper, this is believed (but not guaranteed) to be true. 
However, in the general case these issues must be considered 
carefully. In turn, the goodness of fit enables constraints to 
be placed on simulations of galaxy formation, and can high- 
light missing physics that may be required in order to ex- 
plain the observed properties of galaxies. 



3 MODELS 
3.1 Simulations 

We draw our simulated galaxies from cosmological hydrody- 
nami c simulations run with Gadget-2 (|Springel fc Hernquistl 
120021), including our improvements as described 
in lQppenheimer fc Davi l|2009 . hereafter OD06). This code 
uses an entropy-conservative formulation of smoothed par- 
ticle hydrodynamics (SPH) along with a tree-particle-mesh 
algorithm for handling gravity. Heating is i ncluded via a spa- 
tially uniform photoionizing background (|Haardt fc Madaul 
which is an acceptable approximation for the galaxies 
that are observed at high redshift owing to the fact that 
they form in highly o verdense regions t hat undergo local 
reionization at 2 > 6 (|Dave et al.ll2006ah . All gas particles 
are allowed to cool under the assumption of ionization 
equilibrium, and metal-enriched particles may additionally 
cool via metal lines. Cool gas particles are allowed to de- 
velop a multi-phase interstellar medium via a subresolution 
multi-phase model that tracks condens ation and evapora- 
tion following iMcKee fc Ostrikerl (1 19771 ). Stars are formed 
from cool, dens e gas using a recipe that reproduces the 
iKennicut t (1998) relation; see lSpringel fc Hernquistl (|2003al ) 
for details. The metallicity of star-forming gas particles 
grows in proportion to the SFR under the instantaneous 
recycling approximation. Stars inherit the metallicity of the 
parent gas particle, and from then on cannot be further 
enriched. 

Cosmological hydrodynamic simulations that do 
not include kinetic feedback from star formati o n in- 
variably overproduce star s (e.g., iBalogh et al.l l200ll; 
ISpringel fc Hernquistl l2003al . OD06) . Because superwinds 
can affec t the physical prop erties of the simulated galax- 
ies (e.g. iDave et alj l2006al ). we consider model galaxies 
from simulations with three different superwind schemes: 
(1) a "no wind" model that omits superwind feedback; (2) 
a "constant wind" (cw) model in which all the particles 
entering into superwinds are expelled at 484 km/s out of 
star forming regions and a constant mass loading factor 
(i.e. the ratio of the rate of matter expelled to the SFR) 
of 2 i s assumed (as in the runs of ISpringel fc Hernquistl 
l2003bl ); and (3) the "momentum-driven wind" (vzw) model 
of OD06, in which the imparted velocity is proportional to 
the local velocity dispersion (computed from the potential) 
and the mass loading factor is inversely proportional to the 



velocity dispersion (|Murrav. Quatert. fc Thompson! 20051). 
as inferred from observ ations of local starbursts l|Martinl 
l2005l ; lRupke et al.| [2005'!. This selection is meant to bracket 
plausible models in order to expose any related systematic 
uncertainties; however, owing to the range of successes in 
comparison with IGM metal-line observations obtained 
by OD06 for the vzw model, we focus on this model when 
the conclusions from the different wind models are broadly 
similar. 

All of our wind models were tested in simulations 
that assumed the "old" WMAP-concordant cosmology 
l|Spergel et al.l 120031 ) having n = 0.3, A = 0.7, H = 
70 kms" 1 Mpc" 1 , er 8 = 0.9, and Q b = 0.04. Each of our 
simulations has 2 x 256 3 particles, with parameters as given 
in OD06. We only employ the 16/i -1 Mpc and 32/i _1 Mpc 
simulations from OD06, as the 8/i _1 Mpc runs did not have 
any galaxies large enough to be observable at 2 > 6. An 
additional set of simulations (the "jvzw" model) were run 
using our preferred wind mod el with the 3rd-year WMAP 
cosmology l|Spergel et alll2006j ). namely Q = 0.26, A = 0.74, 
H = 71 kms" 1 Mpc" 1 , cr 8 = 0.75, and Q b = 0.044. Due 
to an error in the initial conditions generation, the power 
spectrum index was set to n = 1 rather than the currently- 
favored n — 0.95; however, this has little impact on our 
results as we will show that they are insensitive to such dif- 
ferences in cosmology. There is a slight change in the wind 
model for jvzw versus vzw, in that jvzw has a smaller mass 
loading factor by a factor of two-thirds compared to vzw 
(in the terminology of OD06, ao = 200 km/s) in order to 
compensate for the lower collapse fraction at high redshift in 
the new cosmology. In addition to 16/i _1 Mpc and 32ft -1 Mpc 
box sizes, we also run a 64/i _1 Mpc box with the jvzw model 
to sample the bright end of the mass function in order to bet- 
ter constrain some observed galaxies that we will consider in 
Sj6] We found that model galaxies from the jvzw simulations 
have bulk properties that are similar to that from vzw. For 
this paper we will compute all luminosity distances assuming 
the new 3rd-year WMAP cosmology. 

We identif y galaxies using S pline Kernel Interpolative 
DENMAX (see lKeres et al.ll2005l for a full description). We 
only consider galaxies with stellar masses exceeding 64 star 
particles, which represents a converged sam ple in terms of 
both stellar mass and star formation history jFinlator et al.1 
2006). According to this criterion, our 16ft _1 Mpc simulation 
volumes resolve galaxies with stellar mass > 1.2 x 10 s M@. 

For this work, the most important output of the sim- 
ulations is the set of SFHs corresponding to the resolved 
galaxies in each simulation at the various redshift outputs. 
We obtain the rest-frame spectrum for each star forma- 
tion event in a given galaxy at the time of observation 
by interpolating to the corre ct metallicity and age within 
the lBruzual fc Chariot! (|2003h models, assuming a Chabrier 
IMF. Summing these up, we obtain the galaxy's intrinsic 
rest frame spectral energy distribution (SED). 

We cons ider the follow i ng pr escriptions for dust red- 
den ing: The Calzetti et al.l (|2000l ) starburst dust screen, 
the ICharlot fc Fall! (120001) embedded star formation law, 
the iGordon et al.l (I2003T) Sm all Magellanic Cloud bar law, 
and the ICardelli et al.1 (119891 ) Milky Way law. We account 
for IGM absorpt ion bluewards of rest-frame L yq using 
the lMadaul l| 19951 ) prescription. The lMadaul (|l995h law may 
be less appropriate for z > 6 than at lower redshifts because 



6 Finlator, Dave, Oppenheimer 



redshift 

30 20 15 12 10 9 B 7 




400 
Time/Myr 

Figure 1. Star formation histories of the 3 best-fitting vzw galax- 
ies for Abell 2218 KESR, along with the 3 best-fitting galaxies 
from each of our one-parameter model SFHs. The SEDs of all 
models match the data with x 2 P er degree of freedom less than 
unity. The vzw SFHs have been sampled in 20-Myr bins and 
smoothed with a 100-Myr tophat for readability. The areas un- 
der the curves are slightly different, reflecting uncertainty in the 
total stellar mass. Note that the best-fitting rising model is very 
similar to the best-fitting simulated model. 



the universe is completin g reionization at this epoch. Indeed, 
ISchaerer fc Pelldl (|2005h found that they were able to im- 
prove the quality of their fits to the SEDs of two reionization- 
epoch ga laxies b y simp ly doubling the optical depth pre- 
dicted by IMadaul ll 19951 ). However, they also found that the 
best-fit derived parameters are relatively insensitiv e to the 
IGM treatment. Thus, for simplicity we retain the IMadaul 
l| 19951 ) treatment without modification. 



3.2 One-Parameter Star Formation Histories 

To date, efforts to use SED-fitting to infer the physical prop- 
erties of high-redshift galaxies have generally employed some 
combination of constant, exponentially decaying, and single- 
burst star formation histories in order to span the presumed 
range of possibilities. In general, it has been found that the 
stellar mass, SFR, and redshift of a galaxy can be fairly well- 
constrained via this technique while the age, metallicity, and 
dust extinction cannot. Much of the gain in precision that 
results from using simulated galaxies in SED-fitting results 
from the relatively small range of SFHs that actually occur 
in the simulations. 

For example, the solid black curves in Figure[T]show the 
SFHs of the 3 galaxies from the vzw simulation that yield 
the best fits to the z ~ 6.7 galaxy Abell 2218 KESR, which 
we will discuss extensively in fJS] The SFHs have been sam- 
pled in 20-Myr bins and smoothed with a 100-Myr tophat in 
order to make the plot more readable. All 3 galaxies begin 



forming stars at z > 15 and exhibit a SFR that is generally 
rising. An examination of simulated SFHs at these redshifts 
shows that steadily rising SFHs are typical. For this rea- 
son we consider a constantly-rising model SFH in this work; 
as we will see, the constantly-rising model reproduces most 
closely the constraints obtained from the simulated galax- 
ies (Figure [8} . This model has to our knowledge not been 
investigated before. 

In order to facilitate comparison with much of the avail- 
able SED-fitting work that is available in the literature, we 
investigate three one-parameter model SFHs for each galaxy 
in addition to simulated SFHs, as described below: 

• Exponentially Decaying SFR We generate models with 
SFR proportional to e _t//T . We use four values of r logarith- 
mically spaced between 10 and 795 Myr, roughly the age of 
the universe for our most distant object. Each of these SFHs 
is sampled at 23 ages t evenly spaced between 10 and 1000 
Myr. 

• Constant SFR We generate models that have been 
forming stars at a constant rate M» for t Myr. For t we sam- 
ple 41 ages that lie between 10 and 1000 Myr, and for M* 
we sample 45 SFRs that lie between 0.2 and 30.0 M Q yr _1 . 

• Constantly Rising SFR In the constantly rising SFH, a 
galaxy's SFR is proportional to its age. While a rising SFH 
can clearly not be maintained for all galaxies until low red- 
shifts, it arises fairly generi cally for high-redshif t galaxies in 
hydrodynamic simulations (|Finlator et al.ll2006T l. We gener- 
ate models in which each galaxy's SFR has been rising at 
a constant rate for t Myr, where for t we have sampled 41 
ages that lie between 10 and 1000 Myr. 

For each star formation history, we have generated mod- 
els with masses in the range log(M,/M0) £ [7.5, 10.5] and 
metallicities Z„/Z Q £ (0.005,0.07,1.0,2.5). These SFHs are 
then put through the Spoc formalism, in order to determine 
the probability distribution of physical properties. During 
the fitting, we require that the oldest star of a given model 
is not older than the age of the universe at the model's red- 
shift. 



4 PERFORMANCE OF Spoc 
4.1 Self-Consistency Test 

We begin by testing that Spoc recovers the (known) proper- 
ties of simulated galaxies. This serves to both test the algo- 
rithm and quantify its intrinsic uncertainties. To do so, we 
take the 73 galaxies that are resolved by our vzw simulation 
at z — 6.5, and determine how accurately we can recover 
their intrinsic physical properties using model galaxies from 
the z — 6 and z = 7 outputs as inputs to Spoc. While the 
model and sample galaxies are not strictly independent in 
this test (all but the least massive galaxies at z — 6.5 cor- 
respond to at least one ancestor in the z = 7 output and 
descendant in the z = 6 output), the galaxies are evolving 
rapidly enough that these populations are effectively inde- 
pendent. The test-case and model galaxies are compared in 
6 bands from i-band to IRAC 4.5/im (the same ones ap- 
plied to Abell 2218 KESR in |5J, where we assume a 0.15 
magnitude uncertainty in each band. The test-case galaxies 
are reddened with a fiducial dust extinction Av = 0.6 via 
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the lCalzetti et al.l (2000) law. We apply Spoc to these test- 
cases using each of the different extinction curves mentioned 
in 33.11 in order to investigate the systematic uncertainties 
resulting from our ignorance of the appropriate extinction 
curve for high-redshift galaxies. During the fitting, redshift 
space is sampled by perturbing each model galaxy over a 
grid extending to Az = 0.5 so that we sample the range 
z G [5.5, 7.5]; Av is sampled over the range Av € [0, 1]. 

SPOC constrains six quantities: M*, SFR, Ay, Z*, age, 
and redshift. The definitions of M* and redshift a r e self - 
evident. Av is defined in terms of the lCalzetti et all (|2000l ) 
reddening presciption. For the purposes of this work, metal- 
licity Z t is defined as the mean mass fraction of metals in 
the galaxy's stars; this is useful in determining what metal- 
licity to choose during population synthesis modeling. Al- 
though metallicity is not the dominant factor in determin- 
ing a galaxy's SED, the fact that the vzw model repro- 
duces the mas s-metallicity rel a tion of star-formin g galax- 
ies at z ~ 2 (|Erb et all feoOfj ; iDave et all 12006b! ') as well 
as for the host galaxy of GRB050904, which is located at 
z = 6.295 (|Berger et al.ll2006l ; iKawai et al.ll2006h . leads us 
to believe that this model's predictions for the metallicities 
of observed reionization-epoch galaxies are plausible (Finla- 
tor et al. 2007, in prep.). We define a galaxy's age as the 
mass-weighted mean age of its star particles; this is more 
meaningful than the more commonly-used age of the oldest 
star, which is both difficult to constrain observationally and 
difficult to predict owing to the stochastic nature of our sim- 
ulations' star formation prescription. We define a galaxy's 
SFR as the average over the last 100 Myr leading up to the 
epoch of observation; if none of a galaxy's stellar mass is 
older than 100 Myr then the age of the oldest star is used. 
This metric is found to correlate more tightly with rest- 
frame UV flux than averages over a shorter time-baseline 
for the numerically simulated SFHs. 

First we consider the case in which the test-case and 
model galaxies are both reddened via the ICalzetti et al.l 
(2000) extinction curve. For this case, the points in Figured 
show how the fractional error in the six inferred proper- 
ties varies with stellar mass and the solid black histogram 
gives their combined distribution. The dotted lines indicate 
the mean formal lcr uncertainties; these are computed di- 
rectly from the probability densities that are returned by 
Spoc rather than from the scatter in the points. In general, 
the recovered physical parameters lie within 50% and 2<r of 
the correct values, suggesting that our SED-fitting technique 
is indeed self-consistent. The fact that the formal uncertain- 
ties are at least as large as the scatter (and, in some cases, 
are somewhat larger) suggests that the formal uncertainties 
are sufficiently conservative. 

The most accurately (and precisely) recovered parame- 
ter is redshift. The high accuracy in this case owes to the fact 
that the /gi4, 2850, and Jno fluxes tightly contrain t he posi- 
tion o f the Lyman break, which itself results from the lMadau] 
(1 19951 ) prescription for IGM absorption. 

Stellar mass is recovered with 20% accuracy, owing pri- 
marily to the fact that the rest-frame optical flux is generally 
dominated by numerous long-lived, low-mass stars whose 
mass-to-light ratio is relatively insensitive to age and dust 
extinction. Additionally, as we will show in Figure [3] the 
lack of a significant systematic offset in the recovered stel- 
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Figure 2. Fractional error in inferred physical properties of 
z = 6.5 vzw galaxies as determined using the z = 6.0 and 
z = 7.0 galaxies as models in SpOC. First we considered the 
case in which both test-ca se and model galaxies are reddened 
via the lCalzetti et al . (2000) extinction curve. For this run, black 
crosses and blue triangles denote test galaxies from the 16/i — x Mpc 
and 32/i — x Mpc volume simulations, respectively; the solid black 
histogram denotes their combined distribution; and dotted lines 
indicate mean fractional la uncertainties computed by SpOC. We 
then compared the same test-case galaxies (reddened with the 
same extinctinction curve) to models with the same SFHs but 
using different dust laws. From these runs, the dotted red, short- 
dashed blue, and long-dashed magenta histograms correspond to 
the cases where the models were reddened with the 



j Cardelli et al.l 

jl989h , iGordon et all ||2003|1 . and ICharlot fc Fall! (120001) laws. 



With few exceptions, the best-fit values are within 50% and 2cr 
of the correct values. Stellar mass and redshift are recovered re- 
markably accurately. 



lar masses owes to the similarity between the SFHs of the 
test-case and model galaxies. 

Metallicity is also accurately recovered. This is ex- 
pected given that there is a tight mass-metallicity rela- 
tion in the simulations (the lcr scatter is 15%) th at does 
not vary strongly with redshift ()Dave et al.l l2006 a,b). and 
the fact that the test galaxies and models came from the 
same simulations. Without this implicit prior, metallicity 
can not be tightly constrained from broadband photome- 
try l|Papovich et alj|200ll ; [Schaerer fc Pelldll2005l ). 

Turning to SFR, we expect a reasonably accurately in- 
ferred SFR given the tight correlation betwee n SFR and stel- 
lar mass that the simu lated galaxies obey l|Finlator et al.l 
120061 ; IDave et al.ll2006al ); in other words, if the redshift is 
known and the stellar mass can be constrained from the 
rest-frame optical flux, then the SFR is already constrained 
to within a factor of two regardless of the rest-frame UV 
flux. Figure[5]bears this out. In detail, SFR is somewhat less 
accurately recovered than stellar mass owing to the degen- 
eracies with age and Av — in fact, a close inspection reveals 
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that galaxies with underestimated SFR have overestimated 
Av and vice-versa. 

Age is accurately recovered owing largely to the small 
range of SFHs that occur in our simulations. Just as only a 
small range of metallicities remains available once the stel- 
lar mass is constrained, a relatively small range of ages is 
available once the redshift and stellar mass are constrained 
(Figured). 

If we relax the assumption that we know the correct 
form of the dust extinction curve, we find systematic ef- 
fects up to the ~ 50% level. In Figure [3j the dotted red, 
short-dashed blue, and long-dashed magenta histograms cor- 
res pond to the cas e s whe re the model s were reddened wit h 
the [Cardelli et all dl989h . SMC bar (|Gordon et alll2003h , 
and Chariot fc F all ( 2000 ) laws while the test-cases were red- 



dened with the Calzetti et al.l l|2000l ) law as before. Metal- 
licity and age are not strongly affected because these are 
tightly constrained by the combination of stellar mass and 
redshift. In contrast, Ay, M», and SFR are underestimated 
for the other curves by up to 60% while the photometric red- 
shifts are systematically off by up to 2%, with the the SMC 
law yielding the largest underestimates. These discrepancies 
owe to the varying slopes of the extinction curves: Steeper 
extinction curves require less overall dust (i.e., lower Av) 
and redder rest-frame UV colors (i.e., lower SFR, and thus 
lower stellar mass in our simulations) in order to match a 
given observed rest-frame UV color. Similarly, photometric 
redshifts are systematically off because the extra suppres- 
sion of rest-frame UV flux that results from an overly steep 
extinction curve can be partially cancelled out by underes- 
timating the galaxy's redshift. 

In summary, stellar mass, metallicity, age, and SFR can 
simultaneously be recovered by Spoc when numerically sim- 
ulated models are used owing to the existence of implicit 
priors on these parameters. Any remaining discrepancy be- 
tween the observed and model UV fluxes is minimized by the 
choice of Av, which is also relatively accurately recovered. 
Thus, our SED-fitting technique is indeed self-consistent. 
However, if the slope of the assumed dust extinction curve 
is incorrect then the resulting best estimates of the physical 
parameters may be off by up to sa 50% while the photo- 
metric redshift may be off by up to 2%. These systematic 
uncertainties are generic to studies of high-redshift galax- 
ies that employ SED-fitting and are unrelated to uncertain- 
ties that result from our ignorance of the correct form of 
high-redshift SFHs. Since it is the latter aspect that we are 
currently trying to constrain, we do not further consider 
SED-fitting errors. 

4.2 Comparison With One-Parameter Models 

It is reassuring but not terribly surprising that Spoc can ac- 
curately recover the physical properties of the galaxies that 
it uses as templates. A more interesting question is how well 
SPOC can recover galaxy properties using a different SFH 
than that of the input galaxy, as this illustrates the varia- 
tions in inferred physical parameters among various assumed 
SFHs. To address this, we have fit the test-case galaxies that 
were used in § 14. II using model sets generated from constant, 
decaying, and rising SFHs as described in § 13.21 

Figure [3] gives the distributions of fractional errors in 
the inferred values of stellar mass, SFR, age, and redshift 
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Figure 3. Comparison of the distribution of fractional errors in 
the inferred physical properties of the galaxies from Figure[2]when 
assuming different SFHs in SpOC. Redshift error is plotted as 
5z/(l + z). All histograms in a given plot have been normalized 
to enclose a constant area. Solid black, dotted cyan, short-dashed 
blue, and long-dashed red curves give the histograms for the simu- 
lated, constant, decaying, and rising model sets; the vertical tick- 
marks at the top indicate the respective medians. Redshift is very 
well-recovered for all models, stellar mass and age to within 50%, 
while SFR can be off by large amounts. 



that result when using the different model sets. The vzw case 
is simply a vertically-binned histogram from Figure [2] Gen- 
erally, the one-parameter models yield stellar mass and age 
results that are within 50% of the correct values. The errors 
for these quantities are generally distributed with slightly 
larger scatter than the errors from the vzw models and show 
systematic discrepancies up to the 40% level. The SFRs are 
overestimated systematically by 50-200% with significantly 
more scatter than returned by the vzw models; this is clearly 
the quantity that is most dependent on the assumed SFH. 
The vzw models systematically underestimate redshift by 
0.2% while the one-parameter models are low by 0.5%; the 
scatters are comparable for all of the models. We briefly dis- 
cuss results specific to each one-parameter model in turn. 

When considering all test-case galaxies together, the 
constant-SFR models tend to underestimate the age and 
stellar mass by 40% and 10% respectively while overestimat- 
ing the SFR by a median factor of 3, the largest discrepancy 
among the SFHs that we consider. When we split the sample 
into "massive" and "low-mass" galaxies at M*/Mq = I0 9 , 
we find that the constant models tend to overestimate the 
ages of "massive" galaxies by ~ 20% while underestimat- 
ing the ages of low-mass galaxies by ~ 40%. In order to 
match the rest-frame optical measurements, the constant- 
SFR models then overestimate the SFR for the low-mass 
and massive galaxies by 100-200% and 0-100%, respectively, 
with 50% scatter in each case. Stellar masses are underesti- 
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mated by 20% for the low-mass galaxies and overestimated 
by roughly the same amount for massive galaxies. This illus- 
trates that uncertainties in parameter recovery are not only 
dependent on the assumed SFH, but also on the mass. 

The decaying models tend to reproduce the stellar mass 
and age with systematic errors of roughly 10% and scatter 
comparable to the scatter from the vzw models. These suc- 
cesses are somewhat surprising because the simulated SFHs 
look nothing like the decaying case. Conversely, the SFRs are 
higher by a median factor of 2.8, only slightly better than the 
constant model. The fact that the SFR could be dramati- 
cally overestimated while the age, stellar mass, and dust red- 
dening (not shown) are recovered accurately probably owes 
to our use of 100-Myr average SFRs. When considering all 
of the stellar mass that has formed in the last 100 Myr, 
a larger fraction of the O-stars will have evolved off of the 
main sequence for decaying or constant SFHs than for rising 
models that show the same 100 Myr average SFR, leading 
to the result that models with differing SFRs nonetheless 
produce similar UV luminosities. The broad success of the 
decaying model despite the input SFHs looking nothing like 
the decaying case shows that SED fitting can yield inter- 
pretations consistent with monolithic collapse models even 
though the true SFH may be quite different. 

The rising models tend to underestimate the age by 
10-40% and overestimate the SFR by 40%, though with sig- 
nificant scatter. Both of these offsets are compensated by 
overestimated Ay in such a way that the stellar masses are 
recovered quite accurately, with < 5% systematic offset and 
scatter comparable to what is achieved via the numerically- 
simulated models. Overall, this model probably recovers the 
true parameters most faithfully among the one-parameter 
models, though it is not a dramatic improvement over the 
others. 

In summary, we have shown that Spoc can self- 
consistently recover the physical properties of the model 
galaxies that we use in fitting observed high-redshift galax- 
ies. Further, simple one-parameter models are able to re- 
cover stellar mass and age to within 50% accuracy and SFR 
to within a factor of three, although there are systematic 
offsets at a comparable level. All models yield photometric 
redshifts with better than 1% accuracy although none of 
them outperform the numerical models. 

4.3 Numerical Resolution 

The Spoc library of simulated galaxies that we employ is 
accumulated from simulations at different volumes and res- 
olutions (e.g. in the ivzw cas e, we use the 16,32,64/i _1 Mpc 
runs ) . In lDave et ail (|2006aT ). we showed that, down to the 
adopted stellar mass resolution limit, the physical properties 
of galaxies are similar at overlapping mass scales between the 
various simulations. To reiterate this point in a way that 
is more relevant to the current work, Figure [3] shows how 
the rest-frame UV-optical color and mean stellar metallic- 
ity vary with stellar mass at z — 6 for resolved (Js 64 star 
particles) galaxies from our 16, 32, and 64 h~ Mpc simula- 
tion volumes. In both cases, the trend and the scatter are 
consistent between the different volumes, giving us further 
confidence that our simulated SFHs are in fact resolved. 

A related way to look for resolution issues is to ask 
whether Spoc will yield similar answers when it is applied to 
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Figure 4. Numerical resolution convergence test I. (Top) Rest- 
frame UV-optical color versus stellar mass at z = 6. The blue 
squares, magenta triangles, and red crosses correspond to the 
resolved galaxies from the 16, 32, and 64 ft _1 Mpc simulations, 
respectively; running medians are also given. (Bottom) Stellar 
metallicity versus stellar mass for the same galaxies. In both cases, 
the median trend and the scatter do not vary with scale, indicat- 
ing that the SFHs of our simulated galaxies do not suffer from 
numerical resolution issues. 



same-mass galaxies at different resolutions. We now demon- 
strate that, indeed, Spoc does recover similar parameters 
for galaxies at different resolutions, and so combining dif- 
ferent resolution simulations into a larger set is justified. Of 
course, one does not need to do so in order to use Spoc, it 
is merely a convenient avenue to increase the dynamic range 
spanned by our model galaxies. 

For convenience, we study A370 HCM6a as it can read- 
ily be fit by our models. We apply Spoc using three sets of 
models derived from the jvzw simulations: once using models 
from the 16/i _1 Mpc volume ( "j 16" ) , once using models from 
the 32/i" 1 Mpc volume ("j32"), and once using both sets. A 
lack of numerical convergence would result in systematic off- 
sets between the probability densities from the first two fits, 
while the combined result shows how the models from the 
two volumes combine to yield our full probability density. 

Figure[5]shows the derived probability density functions 
for the various physical parameters that we consider. The 
ranges agree well, and the best estimates from the j 16 and 
j32 volumes (defined as the means of the probability density 
functions) are consistent at the la level. In detail, the j 16 
models return fits with somewhat lower stellar mass, SFR, 
and metallicity than the j32 models while the j32 models 
yield ages that are younger by about 40 Myr. The small off- 
sets in mass, SFR, and metallicity do not indicate resolution 
problems as they are expected even in the absence of any 
convergence issues. Briefly, the j 16 volume contributes lower 
mass models owing to the slope of the mass function (at the 
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Figure 5. Numerical resolution convergence test II. The dot- 
ted blue and dashed red curves give the probability densities for 
the physical properties of A370 HCM6A as derived from model 
galaxies from a 16h~ 1 Mpc and a 32/i _1 Mpc simulation volume, 
respectively, while the solid black curve results from combining 
the models. All of the probability densities show agreement at the 
lc level and only the age curves show clear evidence of resolution 
effects (see text). This figure indicates that our inferred physical 
properties are not significantly hampered by numerical resolution 



massive end) and our 64 star particle mass resolution cut 
(at the low- mass end). On the other hand, the age offset 
results from a well-known numerical resolution limitation 
whereby galaxies in lower-resolution simulations take longer 
to condense beyond a given critical density in order to begin 
forming stars, yielding younger ages at a given stellar mass 
and redshift. Fortunately, the offset is comparable to the 
intrinsic uncertainty on this parameter. We have repeated 
this test using object SBM03#1 with the 32 and 64/i _1 Mpc 
volumes and found similar results. Hence we do not believe 
that our results using combined simulation samples are sig- 
nificantly hampered by numerical resolution effects. 
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Figure 6. Best-fit spectra from the simulations to the 
data for Abell 2218 KESR, which we have demagnified by 
25 X l lKneib et al.|[20o3) . Spectra from the nw, cw, and vzw, and 
jvzw models are denoted by red dashed, blue dotted, black solid, 
and long-dashed cyan curves, respectively. All four models pro- 
duce galaxies whose spectra match Abell 2218 KESR within the 
errors (\ 2 per degree of freedom < 1). 



veys (ACS) z 85 o, Wide- Field Planetary Camera 2 (WFPC2) 
isi4, and Near-Infrared Camera and Multi-Object Spectro- 
graph (NICMOS) Jiio and Hiqq bands, and only one image 
in the Spitzer /IRAC 3.6 and 4.5 fim bands (the other image 
is blended with a nearby submillime t er so urce at IRAC's 
spatial resolution). Sch aerer fc Pelld |2005) note that the 
fluxes measured by different authors in the optical/near- 
infrared bands disagree due to the inherent difficulties of 
measuring photometry from extended arcs, and that the dif- 
ferent images of the galaxy do not agree in the Hubble/ ACS 
zsso band. Following their suggestion, we use the weighted 
mean of the two images in the optical/near-infrared bands 
(their SED1) and impose a minimum 0.15 mag uncertainty 
in all bands in order to account for differential lensing across 
the images. Q 



5 TEST CASE: ABELL 2218 KESR 



T he triple arc in Abell 2218, dubbed Abell 2218 KESR 
bvlSchaerer fc Pelldl (|2005l ) after its discoverers ijKneib et al.l 
l2004h , is probably the best-studied z > 6 object at present, 
and its physical parameters have been constrained through 
SED fitting by various authors. Hence it provides a good 
test case for exploring the systematics that result from us- 
ing numerically simulated model galaxies, and comparing to 
results employing more traditional simple SFHs. 

The flux from Abell 2218 KESR can be measured from 
two lensed images in the Hubble Advanced Camera for Sur- 



5.1 Modeling Uncertainties: Outflows and Dust 

Figure [5] shows that the SEDs of the best-fitting model 
galaxies from our three galactic outflow recipes and two cos- 
mologies all reproduce the observations with reduced x 2 < 1. 
Moreover, they are remarkably similar. All four models pos- 
sess a very blue rest-frame UV continuum owing to young 
age and low metallicity as well as a pronounced Balmer 

1 ISchaerer &; Pelfdl l l2005h have noted that the published upper 
limits from LRIS and in the Hubble/ ACS Veoe band do not sig- 
nificantly affect the derived parameters. In the case of the Veoe 
limit, we have verified this. 
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SFR (M s yr-') log(M./M e ) redshift 

Figure 7. Probability densities for the physical parameters of 
Abell 2218 KESR based on four numerical simulations assum- 
ing different treatments for superwind feedback. Results from the 
nw, cw, vzw, and jvzw models are given with dashed red, dot- 
ted blue, solid black, and long-dashed cyan lines, respectively. 
So lid green and dotted mag enta boundaries denote r anges inferred 
bv lEgami et al.l [|2005h and lSchaerer fc Pellol [|2005h. respectively 
the l ower limit to stellar mass obtained by ISchaerer fc Pellol 
20051 ) is 7.7). Each curve has been normalized to unit area and 
then scaled so that the three models fit on the same plot. All of 
the best-fit physical parameters except metallicity (see text) are 
robust to changes in the superwind feedback treatment. 

break owing to the presence of older stars. The best-fit pa- 
rameters are similar, indicating that the simulation's ability 
to reproduce the observations is robust to the choice of su- 
perwind feedback prescription and detailed cosmological pa- 
rameters (to the extent of the variations considered). This 
result is akin to the findings from studies employing one- 
parameter model SFHs that good fits can be obtained via a 
variety of assumed SFHs and metallicities. 

To quantify this point, Figure [7] and Table [T] show how 
the derived probability densities for the physical parameters 
of Abell 2218 KESR depend on the choice of wind model. 
The entries in Table [1] list the mean and variance of the 
histograms shown in Figure [7J Each curve in Figure [7J has 
been normalized to unit area and then scaled so that all four 
curves fit on the same plot. 

All of the derived physical parameters except Z* are 
remarkably robust to our choice of superwind feedback pre- 
scription. Generally, the data seem consistent with negli- 
gible dust reddening, a mean stellar age of 100-200 Myr, 
SFRs of 1-3 M Q yr" 1 , a stellar mass of 3-8 x 10 s M©, and 
a redshift z ~ 6.7, in good agreement with other deter- 
minations (|Egami et al.ll2005l ; ISchaerer fc Pe llo 20050- The 

2 Note that because our simulations' star formation treatment 
assumes instantaneous recyling, we give here the total mas s of 
stars that have formed. Using the iBruzual fc Chariot! (|2003T l ta- 
bles, we find that Ri 70% of the stellar mass of stars formed in 



tightness of the constraints on the various intrinsic physical 
parameters results from the relatively narrow range of SFHs 
experienced by the galaxies in the simulations. By contrast, 
the uncertainty in the inferred redshift is determined by our 
treatment of IGM extinction since the inferred redshift is 
dominated by the position of the Lyman-a break. 

Comparing the results from the models in detail, both 
the cw and nw models are less efficient at suppressing star 
formation via outflows. For this reason, at fixed number den- 
sity (or equivalently, dark matter halo mass) the cw and nw 
model galaxies have formed more stars, retain ed more of 
their m etals, have a higher SFR, and are older l|Dave et al.l 
I2006al ). Similarly, at a given stellar mass, the vzw galaxies 
are younger, have expelled a larger fraction of their heavy 
elements, and for these reasons require more dust redden- 
ing in order to match a given observed colour owing to the 
age-met allicity-dust degeneracy. The jvzw models give re- 
sults that are the similar to the vzw results owing to the 
similar feedback treatment. However, the lower values used 
in the jvzw simulation for the cosmological parameters Q m 
and as delay the growth of structure at early times, yielding 
fewer model galaxies for us to match against observations 
for a given simulated volume. As a result, parameter space 
is less well-sampled and the probability density curves (no- 
tably Ay) exhibit more stochasticity. 

Next to stellar mass and redshift, the mean stellar 
metallicity is the most tightly constrained parameter, fol- 
lowed by the mean age and SFR. Metallicity is the only pa- 
rameter for which the different wind models disagree at the 
> lcr level. That the models could agree on all of the derived 
parameters except for the metallicity reflects the fact that 
the effect of metallicity on t he photometry is small c ompared 
to the effect of stellar age l|Schaerer fc Pellol [20051 ) . On the 
other hand, the tightness of the metallicity constraint follows 
from the tight mass-metallicity relation that the simulation 
predicts (§S|. Thus, the apparent disagreement between the 
metallicity constraints is simply a reflection of the tight pri- 
ors imposed by the different simulations combined with the 
robust constraints on stellar mass. Note that because the 
vzw model agrees the best with the distribution of metals 
in the IGM (OD06) a nd the mass-metall icity relationship of 
star- forming galaxies dDave et al.ll2006bJ ). it probably makes 
the most believable prediction of Abell 2218 KESR's mean 
stellar metallicity. 

We investigated the effect that varying the dust pre- 
scription has on t he derived physica l param eters and found, 
in agreement with lSchaerer fc Pellol (|2005l ). that this has no 
significant effect on the derived physical parameters other 
than Ay when the simulated models w ere used. The derived 
Ay were roughly 0.1 mag l ower for the lCardelli et all (|l989l ) 
and I Chariot fc Falll l|2000t ) laws; these differences are ex- 
pected given the extra extinction imposed at rest-frame UV 
wavelengths in the former case and (similarly) for younger 
stars in the latter case. The total amount of light removed by 
dust extinction, which can be regarded as a prediction of the 
total infrared luminosity, is roughly 10 10 Lq, independent of 
the assumed dust prescription. 

Two groups have previously published constraints on 



galaxies that are more massive than 10 Mq at z = 6.75 remains 
in stellar form at z = 6.75. 
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Figure 8. Best-fit combinations of age, extinction, and and SFR 
for 1000 Monte Carlo re-observations of Abell 2218 KESR. The 
solid black, dotted blue, short-dashed red, and long-dashed cyan 
contours enclose 68%, 95%, and 99% of the best-fit solutions for 
the vzw, decaying, rising, and constant model sets, respectively. 
The vzw models allow for the tightest constraints to be placed. 
Interestingly, the 1<t best-fit intervals from the constantly-rising 
SFH are very close to the simulation result. 



this object's properties. lEeami et all (|2005l ) employed a uni- 
form sampling of single-parameter SFHs; their results are 
given in Table [1] and by the solid green vertical lines in Fig- 
ure [7] where we have converted their derived SFR and M* 
to values appropriate for a Chabrier IMF. It is clear that 
their constraints are entirely consistent with our own, al- 
though the tight intrinsic correlations between physical pa- 
rameters in the simulated models allow us to impose tighter 
constraints on all of the derived parameters. In the sec- 
ond work, [Schaerer fc Pellol (120051 ) exhaustively examined 
the systematics of assumptions regarding SFH and template 
spectra with the goal of bracketing the most likely parameter 
space. Their constraints, given by the dotted magenta verti- 
cal lines and listed in[T] are also consistent with ours. Hence 
our simulations, despite a different generic form for their 
galaxies' SFHs than has been assumed in previous investi- 
gations, can reproduce the properties of Abell 2218 KESR 
equivalently well. 



5.2 Comparison to One-parameter SFHs 

Efforts to constrain the physical properties of high- 
redshift galaxies via SED-fitting invariably encounter a 
host of degeneracies between the best-fit parameters, the 
most difficult of which is certainly the age-extinction- 
SFR degeneracy dShaplev et all l200ll ; iPapovich et~aTl l200ll ; 
IShaplev et all l2005h . Generally, young stellar population 
age, high SFR, and low dust extinction all yield bluer photo- 
metric colours while old age, low SFR, and high dust extinc- 



tion all yield redder photometric colours. These degenera- 
cies also contribute to the uncertainty in the inferred stellar 
mass via the stellar age-stellar mass degeneracy, whereby 
older populations have a higher mass-to-light ratio, yielding 
a higher stellar mass at a given measured flux density. Com- 
pounding the problem, there are a number of ways in which 
the best-fit physical properties depend on any assumptions 
that are mad e regarding the shape of a galaxy's star forma- 
tion history (|Shaplev et alj 12005). Numerical simulations, 
in contrast, provide tight relationships between galaxy star 
formation rate, stellar mass, and metallicity (from which 
extinction may be inferred). These predicted relationships 
translate into tighter constraints on physical parameters 
that depend strongly on the shapes of the trial SFHs. 

In order to quantify how effectively the simulations 
reduce well-known degeneracies in best-fit parameters, we 
generated 1000 Monte Carlo re-observations of Abell 2218 
KESR by adding scatter to the photometric measurements 
in a way that was consistent with the reported observational 
uncertainties. For each data set, we then determined the 
best-fit parameters from the vzw models as well as from the 
decaying, rising, and constant SFH models. Figure [8] shows 
the locus of best-fit parameters in several famously degener- 
ate projections of parameter space. The solid black, dotted 
blue, short-dashed red, and long-dashed cyan contours en- 
close 68%, 95%, and 99% of the best-fit solutions for the 
vzw, decaying, rising, and constant model sets, respectively. 

The expected degeneracies that result from one- 
parameter SFHs are easy to see in each panel of Figure [8] At 
the young end, each of the one-parameter models includes a 
parameter space corresponding to a galaxy that has formed 
all of its stars in a rapid burst lasting less than 100 Myr. Our 
use of 100-Myr average SFRs guarantees these models a high 
SFR 4Mq yr" 1 ). Because these models are very young, 
they are intrinsically blue and can require high dust redden- 
ing (Av ^ 1 mag) to match the observations. Additionally, 
their young ages guarantee that their observed optical flux 
is dominated by short-lived O and B stars with low mass- 
to-light ratios, leading to low inferred stellar masses. These 
dramatic burst-dominated models have no analogue in the 
simulations. 

The old end is dominated by the constant-SFR mod- 
els, equivalent to the r — > 00 limit of the decaying models. 
These models are characterized by intrinsically red colours 
and high mass-to-light ratios, leading to low dust extinctions 
and high stellar masses. The oldest of these fits requires the 
galaxy to have been forming stars at ~ lMg yr~ x when 
the universe was less than 10 Myr old; our simulations can- 
not produce this because gas densities have not grown high 
enough to support such SFRs at such early times. 

At the la confidence level, the one-parameter model 
that most closely approximates the simulated galaxies is 
the rising SFH model, albeit with a preference for some- 
what older ages. Returning to Figure [TJ we see that this is 
expected because the simulated galaxies' SFHs are gener- 
ically characterized by a slowly rising SFR at these red- 
shifts. Conversely, at the > 2a level, the constantly rising 
SFH model allows for a wider range of models that have no 
analogues in the simulations. The relatively small range of 
simulated galaxy SFHs typified by the examples in Figure [T] 
leads directly to the relatively tight range of inferred ages 
for the solid contours in Figure[8] Assuming that Abell 2218 
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KESR is located at z = 6.75, the SFHs in Figure Q] suggest 
that it may have formed its first stars before z — 15, with 
(10%, 50%) of its stars in place by z =(11,8); in other words, 
roughly half of the best-fitting models' stars at the epoch of 
observation are over 150 Myr old. This relatively old pop- 
ulation readily accounts for the pronounced Balmer break 
that is visible in Figure [6] and is typical in our simulations. 

In summary, a key point of this paper is demonstrated 
by the fact that, in each panel of Figure [5] the confidence in- 
tervals obtained from simulated galaxies fall well within the 
range described by the complete set of one-parameter model 
SFHs, while yielding the tightest constraints. The tighter 
constraints owe to the relatively small range of SFHs and 
the tight correlations between the various parameters that 
generically occur in hierarchical simulations of galaxy forma- 
tion. This illustrates that, if the simulations are broadly cor- 
rect, the physical properties of high-z galaxies can be more 
precisely constrained using Spoc. However, since at present 
the simulations are larg ely untested at the se epochs (modulo 
the broad successes in iDave et al.ll2006al ). it may be more 
appropriate to regard the tight simulation constraints as pre- 
dictions to be tested against future observations. 

5.3 The Importance of Rest-Frame Optical Data 

In this work we focus on reionization-epoch objects for which 
rest-frame optical measurements are available because the 
rest-frame optical flux is dominated by relatively low-mass 
stars whose mass-to-light ratio evolves slowly relative to 
the stars that dominate the rest-frame UV. For reasonable 
SFHs, this allows tighter constraints to be placed not only 
on the galaxy's stellar mass, but also on its SFH. To quan- 
tify this point, Figure [9] plots the light-weighted median age 
of A2218 KESR versus rest-frame wavelength for the SFHs 
in Figure [1] For example, a point at 2000 A and 15 Myr 
indicates that, for that model, 50% of the photons with rest- 
frame wavelength of 2000 A are emitted by stars that are 15 
Myr old or younger. 

For all of the models that we consider, photons from 
bluewards of the Balmer break are generated by stars that 
are under 100 My r old. This is expected since B stars live 
roug hly 100 Myr (|lbenlll967h and justifies the use of rest- 
frame UV flux as a constraint on the current SFR. Con- 
versely, it explains why rest-frame UV data do not constrain 
a galaxy's SFH prior to w 100 Myr before the epoch of ob- 
servation. By contrast, data from longer wavelengths sample 
the SFH at earlier epochs because the stars that dominate 
these wavelengths live much longer. In the case of A2218 
KESR, the constant and decaying models in Figure [9] sug- 
gest that rest-frame optical data constrain the galaxy's SFH 
roughly 80 and 200 Myr before the epoch of observation, 
while the rising and vzw models fall in between these limits; 
the differences probably owe to the detailed interplay be- 
tween the slope of the SFH and the rate at which low-mass 
stars fade. It is possible that the rough agreement in the op- 
tical portion of Figure [9] explains the tendency of all of the 
best-fit models in Figure [T] to be in rough agreement with 
each other during the interval z m 8 — > 6.7 even though they 
diverge prior to z = 8. Note that this agreement is nontrivial 
given that all of the one-parameter models can in principle 
match the observations with burst-like solutions (Figure [8}. 
Most importantly, it highlights the potential of rest-frame 
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Figure 9. Light-weighted median age of A2218 KESR versus 
rest-frame wavelength. The single solid curve corresponds to the 
best-fitting vzw model while the dotted green, short-dashed red, 
and long-dashed blue curves correspond to the three best-fitting 
constant, rising, and decaying models, respectively. The ranges at 
the top of the figure indicate the full width at 20% of maximum re- 
sponse for an object located at z = 6. The specific one-parameter 
models whose spectra are considered in Figure [9] are the same 
as the models in Figure ITl except that we include only the best- 
fitting vzw model for clarity. Generally, measurements shortward 
of the Balmer break sample light from stars that are 50-100 Myr 
old whereas measurements in rest-frame optical wavelengths con- 
strain older populations. 

optical measurements redwards of the Balmer break to con- 
strain the SFH for reionization-epoch galaxies. 



6 A SAMPLE OF REIONIZATION-EPOCH 
GALAXIES 

While the example of Abell 2218 KESR is illustrative, the 
ultimate purpose of Spoc is to constrain physical parame- 
ters for a large sample of galaxies, in order to characterize 
the observed galaxy population and to constrain the under- 
lying galaxy formation model. To illustrate the sort of in- 
sights gained using Spoc, we apply it to a sample of observed 
z > 5.5 objects that have published broadband photometry 
in the rest-frame UV and optical bands. For each object, 
we used Spoc to determine its best-fit physical parameters 
with the jvzw, constant, rising, and exponentially-decaying 
models. For objects whose redshift has been measured spec- 
troscopically, we run Spoc twice: once with the redshift con- 
strained to the measured value and once with the redshift 
left as a free parameter. The first run allows more accurate 
derivations of the physical parameters, while the second en- 
ables us to study the accuracy of Spoc's photometric red- 
shifts. 

Table [2] compares the means and 95% confidence inter- 
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Figure 10. Spectra of the best-fitting model galaxies overplot- 
ted on the measured photometry. In cases where spectroscopic 
redshifts are available this has been forced; in other cases (Abell 
2218 KESR and Y05#7ab) we used the photometric redshifts. 
Horizontal error bars indicate the full width of each filter at 20% 
of maximum response. The bottom and top horizontal axes give 
the observed and rest-fra me wavelengths. Using the numerically- 
simulated SFHs with the iBruzual fc Chariot! (120031) models, we 
obtain satisfactory fits in all cases and excellent fits in 5 cases. 



Figure 11. Contours enclosing 68 and 99% of the solutions found 
for the various galaxies given the various SFHs that we inspected. 
Thick solid black, dotted red, solid cyan, and dashed blue curves 
correspond to the jvzw, rising, constant, and decaying SFH mod- 
els, respectively. In all cases but GLARE#3001, the jvzw models 
yield results that are generally consistent with the one-parameter 
model SFHs while allowing for equally-tight or tighter uncertain- 
ties. 



vals for the various physical parameters of each galaxy that 
result from the various models, while Figure [11] shows the 
68 and 99% contours in the M*-SFR space for the same 
models. While fitting with one-parameter models, we con- 
fir med that the probabili ty density varies only weakly with 
l|Papovich et al.ll200lD ; for this reason we have omitted 
Z* from Table [2] Our simulations, combined with Spoc, 
suggest that all o f the objects have 0.1 < Z„/Zq < 0.3 
(|Dave et alj|2006af ). 

Figure [TO] compares the synthetic spectra of the best-fit 
jvzw models with the published rest-frame UV through op- 
tical data for our sample. Table[2]shows that we obtain very 
good fits (xt < 2) with the simulated library in all cases ex- 
cept for the unusual object GLARE#3001. Note that these 
spectra and xt values are purely representative; in practice 
we obtain physical constraints through the Bayesian anal- 
ysis described in [JH rather than by simply examining the 
best-fitting model. It is encouraging that the fits obtained 
from our simulated galaxy library are generally equally as 
compelling, in terms of goodness-of-fit measure xtj as the 
fits that result from the one-parameter SFH model libraries. 

Taken together, these findings strongly support the idea 
that most observed reionization-epoch galaxies possess the- 
oretical analogues within our simulations. However, some of 
these objects show minor discrepancies that may be hint- 
ing at some failing in the models. Given the currently large 
observational uncertainties it is difficult to place robust con- 
straints on the models, but it is nevertheless worthwhile to 



discuss each of these systems in detail in order to illustrate 
how Spoc can be used to both reveal physical characteristics 
as well as test the underlying model. 

A 370 HCM 6 A: We have taken the photometric data and 
errors in the observed R,Z,J,H,K' , and S pitzer /IRAC bands 
from ICharv. Stern fc Eisenhardl (|2005l , hereafter CSE05), 
and we have adopted their lensing magnification factor of 
4.5. The measured flux is anomalously large; CSE05 

interpret the excess over the inferred stellar continuum as 
Ha emission line flux and derive an SFR of 140±90Mq yr _1 . 
In our simulations galaxies at z = 6.5 that possess this ob- 
ject's 3.6/^m flux are forming stars at 2-12Mq yr _1 . Since 
such systems are not expected to produce Ha equivalent 
widths that would contribute significantly to the broad- 
band flux, we adopt the uncorrected 4.5/Ltm measurement. 
We verified that, if we use the Ha-corrected 4.5/im flux ob- 
tained by CSE05, the inferred SFR drops from 6.1 ± 2.3 to 
4.9 ± 1.9Mq yr _1 ; this value is in conflict with the assumed 
line strength as expected. The 8.0/im limits quoted in the 
text of CSE05 are different from the limit plotted in their 
Figure 2; we adopt the less restrictive limit from the Figure. 

When we do not enforce the spectroscopic redshift, the 
photometric redshift 5.9 ± 0.4 is dominated by the R — Z' 
limit and is roughly 1.5er below the spectroscopic redshift; 
the formal redshift uncertainty is dominated by the size of 
the R ~ Z' baseline. Although the one-parameter models 
yield more accurate photometric redshifts in this case, the 
difference is comparable to the formal uncertainty and we 
regard it more as a coincidence than as a clue to the SFH 
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of this object. A more accurate photometric redshift would 
require us to include the unpublished measurement in / and 
an accurate Z' transmission curve; we have used the SDSS z 
profile b ecause a transm ission curve for the Z' measurement 
used bv lHu et all (2002) is not available. 

The fits that we obtain are very good, although 
we confirm that the 4.5/^m measurement is difficult to 
interpret as purely stellar continuum emission. In fact, 
even without the anomalous 4.5^im flux the UV con- 
tinuu m is difficult to un derstand owing to the dip at 
K' (|Schaerer fc Pelldl2005l ). and the measured Lya flux is 
difficult to reco ncile with the observed Lya luminosity func- 
tion at z ~ 6.5 (|Malhotra fc Rhoadsi r2004; Kas hikawa et al.l 
2006). Thus, until more precise measurements of this object 
are available, any inferences regarding its physical properties 
should probably be regarded as suggestive. 

With these caveats in mind, we find that using the sim- 
ulated models yields inferred stellar mass, SFR, and Ay 
that are fully consistent with previous derivations as well 
as with the one-parameter model SFHs. The range of stellar 
masses allowed by the jvzw models nicely brackets the range 
between the emission line-dominated and stellar continuum- 
dominated solutions suggested by CSE05, and the r a nge o f 
SFRs is consistent with the derivation of iHu et all (|2002T ) 
based on the Lya line and well below the inferred Ha-based 
SFR of CSE05. Turning to the one-parameter models, the 
widest uncertainties for all of the parameters are provided 
by the decaying model, with the results from the constant 
and rising models spanning a somewhat smaller subset of 
this space. The results from the rising model are generally 
the closest match to the jvzw models although they permit 
a somewhat larger range of SFR; these are burst-like models 
that have no analogue in our simulations. 
SBM03#1, SBM03#3, GLARE#3001: We adopted the 
measured broadband fluxes and spectr oscopic reds h ifts fo r 
these objects from Tables 1 and 2 in lEvles et"aD (2 005). 
Objec t SBM03#1 is the same as object #lab in Yan et all 
(2005). Comparing the derived fluxes in these two papers, 
we find disagreement at the « la level in the ACS 1775, zsso, 
and NICMOS -H160 bands. To be conservative, we therefore 
impose a minimum uncertainty of 0.15 mag in all bands 
for these objects. When fitting for the photometric redshifts 
we includ e the ^775 fluxes; o therwise we exclude this datum 
following lEvles et all (|2005l ). Fi gure [10] shows that we ob- 
tain an excellent fit for SBM03#1 and satisfactory fits for 
SBM03#3 and GLARE#3001, with the chief difficulty in 
th e latter objects bein g the anomalous fluxes in Ks as noted 
by lEvles et all (|2005D . 

For objects S BM03 #1 dDickinson et al.l |2004 
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20031 ) and SBM03#3 jBunker et all 



20031 ; IStanwav et al.N2003T ) Spoc deduces stellar masses of 



~ 10 Mq and mean ages that include the range 100-300 
Myr irresp ective of the assum ed SFH, consistent with the 
findings of lEvles et all ([20051 ). The minimum age for each 
model is older than the ages for the other objects in our 
sample because the small measured uncertainties on the 
3.6/im fluxes for SBM03#1 and #3 create the strongest 
case for a pronounced Balmer break. In our simulations, 
such objects are roughly 130-200 Myr old and are forming 
stars at a healthy 30-60 Mq yr _1 (both la). The resulting 
intrinsically blue colours cause Spoc to select moderate 
dust extinctions Av of 0.4-0.9 and 0-0.1 for #1 and #3, 



respectively, in order to match their UV continuum slopes. 
While these SFRs are within the full range inferred via 
the one-parameter models, they are generally more active 
and younger than inferred via the decaying and constant 
models. The discrepancy owes primarily to the fact that 
the constant and decaying models permit higher SFRs at 
early times (z > 10) than occur in our simulations. Because 
the rising model excludes such early episodes, it yields 
constraints that are more similar to the results from the 
jvzw models. 

The photometric redshifts for SBM03#1 and #3 are 1- 
3a below the spectroscopic values for both simulation and 
one-parameter models, suggesting that unknown systematic 
uncertainties such as uncertainty owing to inaccurate filter 
profiles should be folded into the formal uncertainties; in 
future work an enforced minimum redshift uncertainty of 
5z/(l + z) 0.02 seems reasonable. 

Object GLARE#3001 is a particularly interesting case 
because it represents the worst fit for all our models. This 
object shows a relatively flat SED with little evidence for 
a Balmer break (Figure [TO)) . The photometric redshifts are 
systematically high although they are accurate at the la 
level. Spoc finds that it is perhaps an order of magnitude 
less massive than SBM03#1 and #3. In our simulations, the 
analogues to GLARE#3001 are of roughly the same age as 
the SBM objects while their SFRs are roughly one tenth as 
large, leading to similarly strong Balmer breaks and dust ex- 
tinctions. Comparing with the one-parameter model SFHs, 
Table[2]shows that, for this object, constant and rising mod- 
els yield masses, SFRs, and ages that are marginally consis- 
tent with the results from our simulations (although they 
are skewed to higher SFR) while decaying models yield ex- 
tremely young, burst-like fits whose parameters are entirely 
inconsistent with the jvzw results (this is especially clear in 
Figure [TTJ. 

The chief difficulty in fitting GLARE#3001 is the 
anomalously high flux measured in Ks, as noted by 
lEvles et all l|2005l ). If real, this would suggest that this is 
a low-mass objects undergoing a burst. Table [2] shows that 
the xt P er degree of freedom is lower for the one-parameter 
models than for the simulated models, because the one- 
parameter models have the freedom to produce a higher Ks 
flux through a higher SFR and Ay together with a lower 
stellar mass; galaxies with these combinations of high SFR 
(> IOMq yr _1 ) and low stellar mass (< 10 9 Af©) simply 
do not occur in our simulations. It would be preliminary 
to stake the final interpretation of this galaxy on one broad- 
band measurement (particularly one in Ks), especially since 
even the one-parameter model fits are not particularly good 
(xi > 3). But this object does illustrate how Spoc can pick 
out galaxies that may provide the most stringent tests of 
galaxy formation models. 

Y05#5abc,#7ab: We adopted the Hubble / ACS+NICMOS 
and Spitzer /IRAC fluxe s for object s #5ab c and #7ab from 
the z f» 6 sample in lYan et alJ <|2005l ) and imposed a 
minimum uncertainty of 0.15 mag in all bands as before. 
Spoc failed to find an acceptable fit for object #5abc, for 
either simulated or one-parameter model SFHs. This object 
is evidently a blend of at least three components located at 
different reds hifts and should therefore be fit with multiple 
components l|Yan et al.ll2005l ); however, this is beyond the 
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scope of the present work. We therefore do not show the 
results from #5abc. 

By contrast, we obtain excellent fits for object #7ab. 
This is somewhat surprising given that #7ab is clearly a 
mixture of two components (Figure 1 in lYan et all (|2005h ). 
and probably owes largely to the fact that it does not show 
anomalous single-band fluxes in the way that HCM 6A, 
SBM03#3, and GLARE#3001 do. The photometric redshift 
is constrained to 6.2 ±0 .2, 1.5a higher th an the photometric 
redshift determined bv lYan et al.l (120051 ) (although they do 
not quote an uncertainty). Stellar mass is relatively poorly 
constrained owing to the relatively low signal-to-noise in the 
observed 3.6/im band. As Figure ITTI shows, when we apply 
SPOC with our simulated models, the poorly-constrained 
stellar mass leads directly to a poorly-constrained SFR. 
Comparing with the one-parameter models, we find that 
these span an even larger (overlapping) space. This ob- 
ject illustrates the importance of high-quality near-infrared 
data in order to understand the physical properties of high- 
reds hift galaxie s. 

lYan et all l|2005h applied exponentially-decaying mod- 
els to this object and found a burst-like best-fit solution 
with a stellar mass of log(M*/Mo) = 9.5 (when con- 
verted to a Chabrier IMF), little current star formation 
(~ O.OOlAf© yr~ ), no dust extinction, and an age of 50- 
100 Myr. The stellar mass is fully consistent with our jvzw 
results. However, we find disagreement in the inferred dust 
extinction, SFR, and age because in the simulations galaxies 
of this stellar mass and redshift are invariably older and are 
still forming stars. In particular, we expect Av = 0.34±0.24, 
M» = 9 ± 5M e yr~\ and age = 175 ± 30 Myr (1 a uncer- 
tainty). Once again, the absence of burst-like models in our 
simulations constitutes an effective prior that excludes such 
solutions. When we apply our own exponentially-decaying 
models to this object we find that the probability density 
possesses a sharp peak at low ages (< 100 Myr) and a broad 
but slightly lower plateau at older ages. So while we for- 
mally obtain the lowest % 2 f° r burst-like models, folding in 
our simulation prior leads to an older age being preferred. 
This object is therefore a classic case where simply taking 
the lowest \ 2 among all SFHs results in a substantially dif- 
ferent interpretation than that obtained using physically- 
motivated priors. It will be interesting to see with improved 
obervations whether or not the simulation prediction ends 
up being correct. 

The sixth column of Table [2] combines M» and M* with 
the age of the universe tn at each galaxy's redshift into 
an estimate of the ratio of the present- to past-averaged 
SFR, M,/(M»/tff) = t H /tsFR- A value of t H /tsFR less 
than unity indicates a declining SF H whereas a value greater 
than unity indicates a rising SFHFI lEvles et al. I i|2005l . l2006l ) 
have previously used this quantity to argue that galaxies 
that exhibit Balmer breaks at z ~ 6 must be experiencing 
a declining SFH. In contrast, Table [5] indicates that, while 
this possibility is permitted by the decaying models at the 
2(j level for all but one of the galaxies that we consider, it is 



by no means requirecjf]. For our physically-motivated mod- 
els it is not permitted owing to the generally rising SFHs in 
our simulations, but in fact even for the decaying models it 
is not the favored solution when the full range of possibil- 
ities is taken into account. Hence we do not find that the 
presence of a Balmer break requires a declinin g SFH. The 
disagr eement between our results and those of lEvles et al.l 
(2005) owes partly to the slightly different model sets used 
in the two analyses, and partly to our employing the full 
probability densities rather than restricting our attention to 
the model with the lowest \ 2 ■ Most importantly, however, 
this disagreement emphasizes the fact that the form of the 
correct SFH is not constrained by the current measurements 
and hence the conclusions are highly model-dependent. 

In these fits, we have not considered the possibility of 
any AGN contribution to the SEDs. From rest-UV photom- 
etry alone it is difficult to rule out this possibility. Fitting 
power laws of the form f v oc v a to their rest-frame UV con- 
tinua, we obtain slopes in the range a € [—0.3, —1.1] with 
a weighted mean of -0.61. These a re not sign i ficant ly dif- 
ferent from the slopes inferred by iFan et all (J2001 ) from 
39 quasars at 3.66 ^ z 4.77, with a mean and stan- 
dard deviation of -0.79 and 0.34, respectively. However, 
with Spitzer /IRAC photometry the Balmer breaks in A2218 
KESR, SBM03#1, SBM03#3, and Y05# 7ab are clearly evi- 
dent. Furthermore, IStanwav et ail l|2004al ) and lBunker et all 
(2003) argue that AGN can be ruled out for three of the ob- 
jects (SBM03#1, SBM03#3, and GLARE#3001) based on 
upper limits to the flux of the N V doublet at 1240A, the 
lack of any X-ray detection, and the relatively small veloc- 
ity width of the observed Lya lines (v F whm < 500 kms" 1 ). 
These arguments also apply to A370 HCM 6A. Additionally, 
CSE05 note that the latter object is similar to high-redshift 
Lya emitters, for which the AGN fraction is constrained 
to be less than 5% (e.g., iGawiser et al.| 120061 ; lOuchi et al 
20051 ; iDawson et ail I2OO4I ; IWang et al.l 1200-t ISantos et al 



2004). This suggests that AGN do not dominate their rest- 
frame optical spectra, and emphasizes the usefulness of 
Spitzer /IRAC measurements for constraining the properties 
of reionization-epoch objects. Hence it seems reasonable to 
model all of these objects with stellar population synthesis 
models alone, as we have done. 

In summary, we are usually but not always able to find 
galaxies within our simulations whose photometry suggests 
that they could be the theoretical counterpart to observed 
reionization-epoch galaxies. Using these models as inputs to 
our SED-fitting engine, Spoc, we obtain physical constraints 
on the properties of the galaxies that are, in five out of six 
cases, consistent with the results from one-parameter model 
SFHs. In the remaining case, the observed SED leads the 
one-parameter models to prefer burst-like scenarios that do 
not occur in our simulations. In all cases the tightness of 
constraints from the simulated models is equal to or better 
than that from one-parameter model SFHs. Albeit small, 
this sample of systems illustrates how using Spoc to com- 
pare among SFHs can provide insights into the physics of 



3 Note that, because tfj/tgpn does not involve the inferred age 
of a galaxy's stars, this exercise is not circular; e.g., results from 
fitting exponentially decaying models do not necessarily require 
that tn/tsFR < 1. 



4 Decaying models permit tu/tsFR < 1 while other models do 
not because decaying models can yield blue colors along with a 
lower specific SFR, than the other models. 
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early galaxy formation, and identify unusual galaxies that 
provide the tightest constraints on such models. 



7 CONCLUSIONS 

In this paper we present a Bayesian SED-fitting engine called 
Spoc, which provides constraints on the physical proper- 
ties of galaxies from photometric data. Spoc takes as in- 
put a galaxy's photometry and a set of model galaxy spec- 
tra, obtained either by assuming an analytic form for the 
galaxy's star formation history or from numerical simula- 
tions of galaxy formation, and outputs probability distribu- 
tions for the physical properties of individual galaxies based 
on those model priors. Here we compare and explore im- 
plications for different models, with an eye towards better 
constraining the physical properties of high-redshift galax- 
ies. 

Because Spoc is intended to test whether predicted 
SFHs match those inferred from data, it provides constraints 
on models of galaxy formation that complement compar- 
isons to bulk properties such as galaxy luminosity functions 
or clustering. A successful model must not only reproduce 
observed bulk statistics, but must also reproduce the colors 
of individual galaxies across all observed bands. Galaxies 
with photometry that cannot be well fit in a given model 
provide insights into model failings, and ultimately into the 
physics that drives galaxy formation. Using Spoc it is possi- 
ble to quantitatively determine how well a particular model 
matches an individual galaxy. Such a methodology is partic- 
ularly useful at epochs when only a small number of galaxies 
are detectable, such as for the earliest galaxies known. 

We show that Spoc accurately recovers the input phys- 
ical parameters of model galaxies when fitted with SFH de- 
rived from the simulations themselves, with typical system- 
atic errors less than the formal la fitting errors. When a 
simulated galaxy is fit with one-parameter star formation 
histories or the incorrect extinction model, then larger devi- 
ations can occur. Since the true SFH and extinction law are 
unknown, these deviations can be regarded as systematic 
errors on the determination of physical parameters. In most 
cases, such systematic errors are less than 50% in stellar 
mass, star formation rate, and Ay, and < 2% in redshift. 

We then apply Spoc to six galaxies at z > 5.5 that have 
published photometry spanning the 4000A break (which is 
necessary in order to obtain meaningful constraints on phys- 
ical properties). We begin with a more in-depth study of 
Abell 2218 KESR (z » 6.7), since at present it is probably 
the best-studied reionization-epoch system, and then apply 
Spoc to five more galaxies in order to investigate a wider 
variety of galaxy properties. Our main conclusions are sum- 
marized as follows: 

• The physical parameters derived for Abell 2218 KESR 
using our library of simulated galaxies are in line with 
previous determinations, and those employing various one- 
parameter SFHs. The formal uncertainties are smaller when 
using simulated galaxies, because the simulations intrinsi- 
cally produce a relatively narrow range of SFHs for galaxies 
at these epochs. The existence of an object with the prop- 
erties of Abell 2218 KESR at these epochs does not pose a 
significant challenge to current models of galaxy formation. 



• All inferred physical parameters for Abell 2218 KESR 
except metallicity are remarkably insensitive to the explored 
choices for superwind feedback model, dust extinction law, 
and cosmology. The fundamental reason for this is that re- 
gardless of such choices, simulations produce a similar form 
for the SFH of early galaxies. Disappointingly, this means 
that photometry alone cannot constrain such modeling un- 
certainties. On the other hand, this means that simulation 
predictions of physical parameters using Spoc are quite ro- 
bust. Discrimination between, say, superwind feedback mod- 
els can be obtained through comparisons with b ulk proper- 
ties su ch as luminosity functions (as shown in iDave et al.l 
12006 j ). 

• Exploring a set of six z > 5.5 galaxies, we find that for 
five of the systems the simulated galaxies and one-parameter 
models produce overlapping probability contours and best- 
fit physical parameters that are formally consistent, and the 
minimum reduced y 2 is similar between all models. This 
shows that simulations are usually but not always able to 
produce galaxies with properties similar to observed z ~ 6 
systems. 

• The spectra of simulated galaxies always show a signifi- 
cant Balmer break, despite the fact that their SFHs are best 
characterized as constantly-rising (not decaying). The me- 
dian stellar age is « 120 — 250 Myr (2a) for all six objects, 
spanning a range in stellar masses from ~ 5 x 10 s Mq < 
M* < 2 x 10 10 Mq. Hence the existence of somewhat older 
stellar populations in these early systems is consistent with 
simulation predictions. 

• The object GLARE#3001 is an outlier that is poorly 
fit by our simulations, though in fairness the fits with one- 
parameter models are not optimal either. The best fits tend 
to favor models with younger stellar ages, lower masses, and 
higher SFRs than any simulated galaxies, perhaps indicat- 
ing that this is a galaxy undergoing a burst. However, as 
discussed in the text, these fits are mainly driven by a high 
K s -band flux. If that data point were confirmed it would 
suggest that the simulations may be missing some physical 
process(es) that governs a small fraction of galaxies at this 
epoch. 

Our hydrodynamic simulations of galaxy formation 
make some clear predictions for the star formation histories 
of galaxies at these early epochs. In particular, they predict 
that the stellar mass, star formation rate, and metallicity 
are al l tightly correlated (|Finlator et al.l 120061 ; IDave et al.l 
l2006al ). Indeed such trends are generically seen in most hy- 
drodynamic simulations of galaxy formation, though not 
necessarily in semi-analytic models. Our preliminary com- 
parisons to z ~ 6 galaxies show that, while the simulations 
are not necessarily statistically favored over other classes of 
SFHs, the available data at least do not rule out these sim- 
ulation predictions. Together with the fact that such sim- 
ulat ions can reproduce key bulk properties of early galax- 
ies dFinlator et alj 120061 ; IDave et all l2006aT ) . this suggests 
that models of early galaxy formation are able to reproduce 
a range of observed properties of galaxies at these epochs. 

In this work we have performed the simplest possible 
spectral synthesis calculations, primarily because our focus 
has been a comparison between different SFH models. In the 
future we pl an to investigate diffe rent populat ion synthesis 
models (e.g.. lMaraston et alj|2006h . IMFs (e.g. lFardal et all 
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2006), and nonstellar contributions to the observed fluxes 
such as emission lines from HII regions. Eventually we will 
apply Spoc to a larger sample of galaxies, which will en- 
able a number of interesting investigations. First, the bulk 
statistics (e.g. stellar mass or star formation rate functions) 
derived using Spoc can be compared against that directly 
produced in simulations. While this is partly a circular com- 
parison because simulated galaxies are being used to infer 
the physical properties, in practice there is no guarantee that 
consistency will be achieved because the simulated SFHs are 
quite generic; hence this should provide a stringent test of 
models of galaxy formation. Second, Spoc can be used to 
identify populations of galaxies that deviate dramatically 
from simulation predictions (such as GLARE#3001), in or- 
der to isolate which physical processes may be missing in 
models. In principle this could quantify the contributions 
from nonstellar emission-dominated sources such as AGN 
and extremely dusty objects. Finally, the redshift evolution 
of the galaxy population provides a strong test of models, 
particularly investigating the dramatic change in the na- 
ture of massive galaxi es that appears to occur at z ~ 1 — 2 
jPapovich et al.ll2005l) . 

SPOC is a general-purpose code, in principle able to uti- 
lize any type of model that produces detailed galaxy SFHs, 
be it a hydrodynamic, semi-analytic, or analytic model. In 
principle there is nothing that restricts Spoc to high red- 
shift use, but in practice other physical phenomena such as 
AGN contamination may become more important at lower 
redshifts. Spoc could also provide a useful tool to identify 
galaxy classes that are not readily reproduced in current 
models. We are currently planning to make Spoc and our 
latest library of simulated galaxies publicly available upon 
publication of this paper. We hope that Spoc will be a use- 
ful and flexible tool for conducting detailed comparisons be- 
tween simulations and observations, as are critical for ad- 
vancing our understanding of galaxy formation. 
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Table 1. Derived Parameters for Abell 2218 KESR 



Model a 


A v b 




Age (Myr) 


MAMq yr" 1 ) 


log(A/,/M Q ) d 


z 


nw 


0.21 ±0.14 


0.0041 ± 0.0007 


147 ± 30 


2.2 ±0.6 


8.71 ±0.13 


6.71 ±0.06 


cw 


0.14 ±0.12 


0.0035 ± 0.0005 


160 ± 33 


2.0 ±0.6 


8.73 ±0.14 


6.71 ±0.06 


vzw 


0.17 ±0.12 


0.0018 ± 0.0002 


175 ± 25 


2.2 ±0.6 


8.78 ±0.12 


6.71 ±0.06 


jvzw 


0.23 ±0.13 


0.0020 ± 0.0002 


160 ± 12 


2.4 ±0.7 


8.74 ±0.12 


6.71 ±0.06 


Egami et al. (2005) 






50-450 


0.1-3.3 


8.5-8.8 


6.6-6.8 


Schaerer & Pello (20051 


^ 0.2 




3-400 


1.5-1.9 


7.7-9.0 


6.0-7.2 



a nw, cw, vzw denote no, constant and momentum-driven winds, respectively; jvzw uses 3rd-year WMAP cosmology. 
b Assumes the lCalzetti et all l|2000l ) dust prescription and Rv =4.05. 
c metal mass fraction in stars 

d Assumes a Chabrier IMF. To convert from a Chabrier IMF to a Salpeter IMF from 0.1-100 M Q , multiply by 1.5. 
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Table 2. Mean and 2a Parameter Ranges for z 5.5 Galaxies 



Model 3 A v b Age/Myr M* C (M yr" 1 ) log(M»/M Q ) c t H /t SF R d z c X V 



A2218 KESR 

jvzw 0.23 (0-0.53) 160 (132-181) 2.4 (1.2-3.7) 8.7 (8.5-9.0) 3.5 (2.5-4.2) 6.71 ± 0.06 0.57 

constant 0.42 (0-1.0) 107 (7-350) 6.7 (2-25) 8.7 (8.1-9.1) 17.2 (1.5-86.1) 6.68 ± 0.08 0.62 

decaying 0.34 (0-1.00) 104 (8-370) 6.9 (1.0-25.2) 8.7 (8.1-9.1) 19.2 (0.5-86.0) 6.68 ± 0.07 0.55 

rising 0.33 (0-0.87) 129 (5-268) 3.5 (1.5-12.5) 8.7 (8.2-9.1) 6.7 (2.2-56.3) 6.69 ± 0.07 0.60 

A370 HCM#6A 

jvzw 0.49 (0-1.08) 191 (122-250) 7.8 (2.7-18.2) 9.3 (8.9-9.6) 2.9 (2.3-4.4) 6.56/5.87 ± 0.35 1.7 

constant 0.76 (0.13-1.24) 182 (24-402) 14.0 (4.1-28.1) 9.6 (9.0-9.9) 3.9 (1.1-17.3) 6.56/6.18 ± 0.41 1.5 

decaying 0.55 (0-1.37) 217 (9-502) 16.4 (0.6-83.5) 9.5 (8.7-9.9) 9.2 (0.4-86.1) 6.56/6.15 ± 0.41 1.5 

rising 0.82 (0.13-1.37) 138 (8-282) 16.1 (4.2-40.1) 9.4 (8.8-9.8) 6.1 (2.2-43.4) 6.56/6.15 ± 0.41 1.6 

SBM03#1 

jvzw 0.67 (0.36-0.93) 175 (123-231) 42.1 (30.9-62.1) 10.1 (9.9-10.2) 3.7 (2.4-4.7) 5.83/5.64 ± 0.08 1.7 

constant 0.45 (0.25-0.62) 388 (235-498) 24.5 (18.0-29.3) 10.2 (10.1-10.4) 1.4 (1.0-2.1) 5.83/5.68 ± 0.08 1.2 

decaying 0.29 (0-0.75) 357 (89-601) 20.7 (0.63-99.6) 10.3 (10.0-10.5) 1.3 (0.1-10.1) 5.83/5.68 ± 0.08 0.9 

rising 0.78 (0.50-1.12) 212 (51-332) 51.7 (29.5-102.0) 10.2 (10.0-10.4) 3.4 (2.0-9.0) 5.83/5.65 ± 0.08 1.6 

SBM03#3 

jvzw 0.44 (0.09-0.69) 170 (124-231) 46.1 (24.9-75.3) 10.1 (9.9-10.3) 3.9 (2.4-4.8) 5.78/5.68 ± 0.07 2.2 

constant 0.18 (0-0.37) 385 (225-498) 25.0 (18.0-29.3) 10.3 (10.0-10.4) 1.4 (1.0-2.3) 5.78/5.74 ± 0.07 1.7 

decaying 0.25 (0-0.87) 301 (30-601) 41.9 (8.3-251.0) 10.3 (9.9-10.6) 3.0 (0.2-20.4) 5.78/5.71 ± 0.08 1.5 

rising 0.58 (0.25-1.00) 188 (25-332) 64.9 (28.4-159.2) 10.2 (9.9-10.4) 4.4 (2.1-13.7) 5.78/5.69 ± 0.08 2.0 

GLARE#3001 

jvzw 0.19 (0-0.46) 183 (123-253) 5.2 (3.4-8.7) 9.2 (9.0-9.3) 3.7 (2.4-5.9) 5.79/5.86 ± 0.07 8.7 

constant 0.47 (0.12-0.62) 31 (6-126) 18.9 (6.2-29.3) 8.8 (8.4-9.2) 37.3 (4.5-101.2) 5.79/5.81 ± 0.07 3.5 

decaying 0.59 (0.13-0.87) 15 (6-85) 30.3 (7.4-63.3) 8.7 (8.4-9.2) 63.5 (5.6-101.2) 5.79/5.80 ± 0.08 3.3 

rising 0.53 (0-0.87) 42 (5-200) 17.6 (5.1-40.1) 8.8 (8.4-9.3) 38.1 (3.4-101.2) 5.79/5.81 ± 0.08 3.0 

Y05#7ab 

jvzw 0.34 (0-0.94) 175 (126-227) 9.2 (3.7-26.7) 9.4 (9.0-9.8) 3.4 (2.6-4.8) 6.25 ±0.16 0.56 

constant 0.56 (0-1.00) 210 (27-448) 16.2 (4.7-28.7) 9.7 (8.8-10.3) 4.3 (1.2-19.0) 6.17 ±0.17 0.18 

decaying 0.48 (0-1.24) 212 (10-532) 22.7 (0.8-105.3) 9.7 (8.5-10.4) 12.4 (0.5-92.8) 6.15 ±0.18 0.14 

rising 0.63 (0-1.37) 140 (5-299) 20.3 (4.6-63.8) 9.5 (8.7-10.2) 7.5 (2.2-61.1) 6.18 ±0.18 0.35 



a "jvzw" denotes the simulated galaxies; "constant", "decaying", and "rising" denote representative one-parameter model SFH models. 
b Assumes the lCalzetti et ail l|2000l ) dust prescription and Ry = 4.05. 

c Assumes a Chabrier IMF. To convert from a Chabrier IMF to a Salpeter IMF from 0.1-100 M Q , multiply by 1.5. 
d The ratio of the current SFR to the past-averaged SFR, defined as tu/tsFR 

c Where 2 numbers are given, the first is the spectroscopic redshift and the second is the photometric redshift. 

f Minimum reduced \ 2 f° r each combination of models and observed object. The number of degrees of freedom is given by the number of 
detections used minus 3 for A2218 KESR and Y05#7ab, and minus 2 for the others. 
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